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Absence of the Kasner singularity in the effective dynamics from loop quantum 

cosmology 

Ghanashyam DatcQ 
The Institute of Mathematical Sciences, CIT Campus, Chennai-600 113, INDIA 

In classical general relativity, the generic approach to the initial singularity is usually understood 
in terms of the BKL scenario. In this scenario, along with the Bianchi IX model, the exact, singular, 
Kasner solution of vacuum Bianchi I model also plays a pivotal role. Using an effective classical 
Hamiltonian obtained from loop quantization of vacuum Bianchi I model, exact solution is obtained 
which is non-singular due to a discreteness parameter. The solution is parameterized in exactly 
the same manner as the usual Kasner solution and reduces to the Kasner solution as discreteness 
parameter is taken to zero. At the effective Hamiltonian level, the avoidance of Kasner singularity 
uses a mechanism distinct from the 'inverse volume' modifications characteristic of loop quantum 
cosmology. 

PACS numbers: 04.60.Pp,98.80. Jk,98.80.Bp 



While the celebrated singularity theorems of classical 
general relativity imply that the backward evolution of an 
expanding universe leads to a singular state, the nature 
of the singularity is elucidated in terms of the Belinskii- 
Khalatnikov-Lifschitz (BKL) scenario Q, The sce- 
nario views the spatial slice close to the (space-like) 
singularity as made up of approximately homogeneous 
patches each of which evolves according to the vacuum 
Bianchi IX model. As the singularity is approached, the 
patches fragment indefinitely, asymptotically becoming 
infinitely small. The vacuum Bianchi IX evolution in 
turn, can be viewed as a succession of Bianchi I evolu- 
tions (Kasner epochs), interleaved by transitions among 
the Kasner epochs. This interleaved evolution continues 
indefinitely in a chaotic manner. Apart from the indef- 
inite number of Kasner epochs and transitions among 
these, the monotonic decrease of the volume and the sin- 
gular nature of the Kasner solution are responsible for 
the infinite fragmentation of the homogeneous patches. 

Singularities 'evolving' from non-singular physical sit- 
uations are indicative of breakdown of the extrapolation 
of the classical evolution i.e. dynamics of Einstein equa- 
tion and call far an extension/modification of the classi- 
cal theory/framework. One natural avenue is to appeal 
to a corresponding quantum theory of gravity. Since the 
classical picture of the singular behaviour involves highly 
dynamical geometries with arbitrarily large curvatures, 
a quantum theory which does not depend on any pre- 
selected background geometry is likely to be most suit- 
able for obtaining the required extensions. 

Loop quantum gravity (LQG) is precisely such a back- 
ground independent approach While the LQG for 
the general inhomogeneous situations is quite compli- 
cated, its methods can be implemented and tested in 
simpler contexts of spatially homogeneous geometries. 
Indeed, loop quantization of the so-called diagonalizcd, 
Bianchi class A models has been carried out (j, [j| and 
shown to be non-singular within the quantum frame- 
work. In the quantum framework, non-singularity means 



non-breakdown of the fundamental dynamical equations 
(which are partial difference equations) and boundedness 
of relevant operators such as the inverse triad operator 
which enter the quantization of curvature components, 
matter densities etc. 

It is of course more intuitive and convenient to obtain 
the modifications implied by the quantum theory in the 
familiar geometrical setting of classical general relativity 
i.e. obtaining the modifications to the Einstein dynamics 
keeping the kinematical framework of Ricmannian geom- 
etry in tact. This has been done systematically for the 
isotropic models 0,0,13 in terms of an effective Hamilto- 
nian. The derivation of the effective Hamiltonian is based 
on the observation that if the fundamental dynamical 
equations admits a solution which is WKB approximable 
(i.e. the amplitude and the phase are slowly varying 
in a suitable sense) in some domain, then to the order 
H° one obtains a Hamilton- Jacobi equation from which 
an effective classical Hamiltonian can be read-off, also 
valid within the same domain. The largest possible do- 
main of validity of such an approximation is constrained 
by the classical 'turning points' dictated by the effective 
Hamiltonian. The form of the effective Hamiltonian so 
obtained does not depend on details of the presumed so- 
lution, but of course the actual domain of validity of the 
WKB approximation and consequently of the effective 
Hamiltonian is sensitive to the presumed solution. Being 
o(H°), the effective Hamiltonian is insensitive to factor 
ordering issues. 

There are various types of corrections that arise. The 
most dramatic one is the correction implied by the non- 
trivial quantization of inverses of various classical quan- 
tities such as scale factors, triad components, volume etc 
|3, llfj ■ In the isotropic context, this has lead to a vari- 
ety of implications UJI llM LL-I ■ For the anisotropic 
context, this leads to suppression of chaotic behaviour of 
the Bianchi IX model [15j with the further result that 
asymptotically for vanishing volumes, a Bianchi IX so- 
lution approaches a Kasner solution in a stable manner 
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is an order 1 parameter and, 



For the vacuum Bianchi I model, there is no potential 
term and no scope for a modification of dynamics due to 
quantizations of inverses of triad components. Although 
at the quantum level, there is no singularity |^|, the ef- 
fective classical Hamiltonian derived from a continuum 
approximation, is identical to the Einstein Hamiltonian 
leading to the same singular Kasner solution. Recently 
however an alternative method of deriving effective clas- 
sical Hamiltonian has been devised, in the context of 
isotropic models |8( , which uses the WKB ansatz directly 
at the difference equation level bypassing the step of first 
deriving the Wheeler-De Witt differential equation to be 
followed by WKB approximation. The same method can 
also be applied in the anisotropic context which leads 
to an effective Hamiltonian different from the Einstein 
Hamiltonian. Exact, non-singular solution of this Hamil- 
tonian for the vacuum Bianchi I model is the result pre- 
sented here. 

A Bianchi I space-time is specified in terms of the met- 
ric of the form, 



ds 2 = alt 2 -^a^idx'f 



(1) 



where t is the synchronous time. The vacuum Einstein 
equations then lead to the well known Kasner solution: 
aj(t) ~ t 2ai where ai are constants satisfying the con- 
ditions Y)j a j = 1 = ^2 j a i- F° r subsequent comparison 
with the new solution, it is convenient to describe the 
time evolution in terms of a new time coordinate r corre- 
sponding to the lapse J\f :— 010203, defined by J\fdr = dt. 
The scale factors then evolve as ai ~ e a ' T . 

Loop quantization of all diagonalized, Bianchi class 
A models has been given in Briefly, it may 

be summarized as follows. The kinematical Hilbert 
space is spanned by orthonormalized vectors labeled as 
|Mi> M2> M3)) i 1 ! £ ^- These are properly normalized 
eigenvectors of the triad operators p 1 with eigenvalues 
^j£pfii, where 7 is the Barbero-Immirzi parameter and 
£p := 8irGh :— kH. The volume operator is also di- 
agonal in these labels with eigenvalues V(fl) given by 
(^7^p) 3 ^ 2 '\/l/ i iM2M3|- Here we have used the vector no- 
tation to denote the triple (^1,^2,^3). Imposing the 
Hamiltonian constraint operator on general vectors of the 
form |s) = s(/2)|/2) leads to the fundamental differ- 
ence equation for the coefficients s(fl). Here the sum is 
over countable subsets of K 3 . There are further gauge 
invariance conditions |J] which do not concern us here. 

In the present context of vacuum Bianchi I model, the 
fundamental difference equation takes the form 

V* A 12 (fl; e 12 )s(fi; 6*12) + cyclic = , where (2) 



e*i2 = (ei, £2, ^2) with each of the e* taking values ±1; 
s(/2;e*i2) = s(/xi - noex - ^[,^2 - Mo£2 ~ ^2,^2); Mo 



-4i2(/2; £12) 



V(fl; ei2)d(^ 3 )(eie 2 + ^2) 







(3) 

/! = 



In the summary above, we have used the non-separable 
kinematical Hilbert space and also made the parameter 
yit explicit 

To derive the effective Hamiltonian, one assumes that 
there exist solution(s) of the partial difference equa- 
tion which have a slowly varying amplitude and phase 
at least in some region of large volume. Explicitly, 
defining p(f-ii) := ^jipfl, one introduces an interpolat- 
ing function, ip{p)i such that s(/2) — ip(p(pi)) = C{p) 
exp{j-$(p)} and assumes that the amplitude and phase 
are slowly varying functions of p in the sense that higher 
order terms in the Taylor series about any p in the rel- 
evant region, are smaller than the lower order terms 
when compared over the quantum geometry scale q := 
Taylor expanding the interpolating wave func- 
tion i/i(p; e/j), (IJ) — (12, 23, 31), it is straightforward to 
check that the leading terms (in powers of h) in the real 
part of the equation are o(h°) while those in the imagi- 
nary part are o(h). The o(fi°) terms involve only the first 
order partial derivatives of the WKB phase. Identifying 
K\ := K^\^t^(j)) etc., one infers the Hamiltonian system 
from the Hamilton- Jacobi equation as: {p 1 ,Kj} = kSj, 
with the Hamiltonian (e :— /io7), 



nAfH cS (p,K) = -2 



-, ? sinei^i sinei^ ,. 
p p 1- cyclic 



(4) 



e will be referred to as the discreteness parameter. 

The effective Hamiltonian is periodic in eKj and we 
may restrict our attention to — 7r < eKj < tt. The 
imaginary part of the equation however requires the 
domain of validity to be restricted further in order to 
be self consistent with the assumption of slow varia- 
tion of the interpolating wave function. The restric- 
tion is: pi > q and small neighbourhoods of ±^ are 
to be excluded. Thus along the eKj axes, the effective 
Hamiltonian is a good approximation in the intervals: 
(-7T,-f - £),(-§ + 8, § - 5),(| + 5,7r) for some small 
positive 8. The values eKi = ±f will turn out to be the 
'turning points' of the trajectories, p(r), of the effective 
dynamics. 

For eKj ~ 0, one can use e _1 sin eKj sa Ki. Then the 
effective Hamiltonian goes over to the Einstein Hamilto- 
nian and the dependence on e drops out. Keeping only 
the H° terms, we have effectively taken \p T \ larger than 
the quantum geometry scale set by q. Thus it is clear that 
Einstein dynamics is reproduced for large values of triads 
and small values of their conjugates Ki. In this regime, 
the Hamilton's equations for p 1 , combined with the rela- 
tion of the triad components to scale factors, \p ! \ — ojok 
(and cyclic permutations), identifies the Kfs as compo- 
nents of the extrinsic curvature of the constant t slices: 
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Kr 



1 da i 
' 2 dt 



The effective Hamiltonian thus deviates 



from the Einstein Hamiltonian mainly in the region of 
the phase space with not too small values of tKj. 

The small value of eKi can be achieved by taking the 
limit e — > [Kj fixed) which then removes any restric- 
tion on Kfs. From a purely classical perspective, such 
a view may be welcome especially since the parameters 
/io,7 disappear which are absent in the Einstein dynam- 
ics. However, from a quantum perspective, the discrete- 
ness parameter cannot be zero |l7j . For a non-zero value 
of discreteness parameter, the effective Hamiltonian rep- 
resents an extension of Einstein dynamics for phase space 
regions beyond small extrinsic curvatures and large triad 
components. Since the non-zero value of e reflects a dis- 
crete structure in a specific technical sense, the effective 
dynamics is to be viewed as a way of extending the Ein- 
stein dynamics by incorporating the discrete feature of 
quantum geometry. 

We use geometrized units (k = 1 = c) so that all quan- 
tities have dimensions of powers of length. The triad 
variables p 1 and h have dimensions of (length) 2 ; the ef- 
fective Hamiltonian, scale factors, the synchronous time 
all have dimensions of length while the r has dimensions 
of (length) _2 . Using the naturally available quantum ge- 
ometry length scale of yfq, all dimensionfull quantities 
below, are rendered dimensionless. We will continue to 
use the same symbols though. 

The Hamilton's equations from the effective Hamilto- 
nian are, 



dp 1 
~dr~ 
dKi 

~dV 
o 



-2p coseKi 

-2 



p J sineKj p K smeKx \ 



e e 
s'meKj f p J smeKj p K s'meKK 



cyclic 



e \ e 
p 1 sineKi\ f p 2 smeK2 



) ,(5) 
(6) 
(7) 



— are con- 
= £/ aj - 



It follows immediately that p' sir ^ Kl 
stants of motion, with {ctiotz -I- cyclic) 

(J2i a i) 2 following from J7J). 

If all the ai are zero, then all the p 1 , Kj are also con- 
stants. Since lapse is non-zero, the p 1 are non-zero con- 
stants and the solution represents the usual Minkowski 
space-time. The effective dynamics thus retains the 
Minkowski space-time as a solution indicating a good 
classical limit of the quantum dynamics. For a non-trivial 
solution, then, J2i a ? ^ must hold and as usual, by a 
constant scaling of the N (or of r) , one can arrange the 
ai to satisfy: Y) j cti = 1. Thus the the one dimensional 
parameter space of these solutions is exactly same as that 
of the usual Kasner solutions. (The special cases of the 
form ati = 1, a.% = 0:3 = are independent of e in their 
behaviour and are not considered here. Thus all ai are 
assumed to be non-zero.) It follows that exactly one of 
the a/'s must be strictly negative and the remaining two 
strictly positive and aj + olk = 1 — a/ > 0. 



Since ai are non-zero constants, neither p 1 nor sineif/ 
can vanish and thus cannot change sign. For dcfiniteness, 
let us restrict our attention to 'positively oriented' triad 
and in particular take p 1 > 0, V7. Then, sgn(sineifj) = 
- sgn(aj), which fixes the two quadrants to which the 
'angles' 8] :— eKj must be confined along a solution 
i.e. either < eKi < ir or — ir < eKj < 0. Clearly, 
p 1 approach 00 as eKi approach the end points of the 
intervals and take the minimum value • for eKj = 
±-|. The traversal of 9i with r is also fixed by JSJh 0[ 
must decrease if sin#j is positive and increase if sin#j is 
negative i.e. we must have Qj — > 0± as r — > 00 (8j — > ±7r 
as r — > —00). 

Eliminating K] in favour of a / , p 1 , leaves us with an 
equation for p 1 , namely, 



dp 1 
~dT 



= ±(l-a 7 W(p') 



\ 2 J 



(8) 



The ± is determined by the quadrant to which the angle 
eKj belongs. The solution is easily obtained as, 

At) = cosh {(l - «i)(t - (rr)o)} • (9) 

Notice that a triad component attains its smallest value, 
^Y^, at t — (17)0 while for large \t\ it behaves as p 1 ~ 
(e|aj|/4)exp((l-aj)H). 

In terms of the scale factors, aj = vo ^ mc | the solution 
is given by, 



a?(r) 



1 — Oil 



(cosh{(l - aj)(T - (rj) )}) X 
(cosh{(l - a K )(r - (t k )o)}) x 
(coshia-a^^-^^D^llO) 



For comparison, the triad and the scale factor for the 
Kasner solution are, 



p I (r)=p^ 1 -^ , a]{r) = {a I )le 



2 2ajr 



(11) 



One can recover the Kasner solution from the modified 
one by taking the limit (tj)o — > —00, e — > such that 



f m\ e -{l-ai)(Ti) 

e 4 e 



For the Kasner solution, two scale factors vanish and 
third one diverges such that the volume vanishes expo- 
nentially with t, as t — ► —00. This translates into a 
finite synchronous time t in the past, making the Kasner 
solution singular. 

By contrast, for the modified solution, none of the 
triad variables can become zero at any r and the vol- 
ume never vanishes. Furthermore, due to the hyperbolic 
cosine function, for both asymptotic times, the modified 
solution approaches the large volume behaviour of the 
Kasner solution. In particular, the scale factor behaves 
as aj — » (e(l — aj)/4)exp(2aj|r|) as |r| — * 00. Conse- 
quently, exactly one scale factor vanishes while the re- 
maining two diverge as |r| — > 00. Thus if one begins with 
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a cubical cell at some finite r, then the cell will become 
planar after both forward and backward evolution (under 
Kasner evolution the cell will become planar in 'future' 
and one dimensional in the 'past)'. Since r — > ±00 cor- 
respond to t — > ±00, the vanishing/diverging behaviour 
of scale factors never occurs for finite t. The modified 
solution is thus non-singular. 

Consider the behaviour of the volume, V 2 — p 1 p 2 p 3 . 
For |r| — > 00, V 2 — > §4 |a!ia!2a3|e 2 ' r ', and therefore it 
must have a minimum, V*, for some r*. It is clear that 
for t larger (smaller) than all (tj)o, V^ 2 is monotonic. 
Thus r* must lie in the interval of the minimum and 
maximum values of the parameters (tj)q. The minimum 
volume, V*, reached by any particular solution depends 
on a as well as on tq. For any given df, the smallest 
possible minimum volume is attained for solutions for 
which all (77)0 are equal. Parameterizing the a as 
a\ := — u/(l+u+u 2 ), a2 := (l+u) / (1+u+u 2 ) , := (u+ 
u 2 ) I (1 + u + u 2 ) , u > 1, this smallest possible minimum 
volume for a given u, is given by, 



(K)r 



|aiQi2Q3| 



u(l + u) 
(1 + u + u 2 )^ 



(12) 



which can be arbitrarily small though strictly positive. 

Since singularity theorem is evaded, 'energy' condi- 
tion^) must be violated. In the present context, this 
means that R u < must hold some where in the space- 

time. By computing R u — ~J2i a 7 ~d§^> one can see 
easily that this is so in the neighbourhood of V* . 

So far we have focussed on the features of the exact so- 
lution. Recall that for the validity of the effective Hamil- 
tonian we also need p 1 > q, or in the dimensionless vari- 
ables used above, p 1 > 1. This puts a restriction on the 
time t for which any specific solution can be trusted as an 
approximation. Specifically, cosh{(l — or)(r — (tj)o)} > 
7j§7j- Had we used the effective Hamiltonian obtained 
from the continuum approximation 0, Hfij | which is same 
as the Einstein Hamiltonian, we would get the usual 
Kasner solution and obtain the restriction on r as, exp 
{(1 — a/) (r — (tj)o)} > elaTJ' ^or ease °^ com P ar i son > 
we have just written pg := ^^-!-exp{— (1 — a/)(r/)o} so 
that at for large t, both the modified and the usual 
Kasner solutions match. Since the modified solution 
is always larger than the Kasner solution at the same 
t, the t— regime of validity for the effective solution is 
larger than that for the Kasner solution. However, since 
> I > 2, the parameters (r/) are out side the do- 
main of validity and so is the smallest (non-zero) value 
of triad components. The 'bounce' in the triad com- 
ponents is thus in the quantum regime. This is differ- 
ent from the isotropic case ^3 where the bounce im- 
plied by effective dynamics can occur in the domain of 
validity of the effective Hamiltonian, depending upon de- 



tails of matter Hamiltonian. Thus, in the simplest of the 
anisotropic models, although the effective dynamics in- 
ferred is non-singular due to bouncing triad components, 
these bounces lie in the quantum domain and for a 'reli- 
able' removal of singularity, one still needs to appeal to 
the quantum theory. 

In summary, we make three points: (1) the method of 
effective Hamiltonian can be extended to homogeneous, 
anisotropic models and leads to non-singular effective 
dynamics with the exact solution parameterized by the 
same Kasner parameters; (2) unlike (at least some of) the 
isotropic models, for reliable singularity removal, one has 
to appeal to quantum theory; and (3) the quantum the- 
ory stipulates modifications of the classical Einstein dy- 
namics, not only for small triads (p 1 ) but also for larger 
values of their conjugates (K]) which is responsible for 
the non-singularity. 

I would like to thank Martin Bojowald and Golam Hos- 
sain for useful remarks. This work was initiated during 
my visit to AEI, Golm during December 2004. The warm 
hospitality is gratefully acknowledged. 
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